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Scale setting for QCD with two flavours of staggered quarks is examined using Wilson flow over 
a factor of four change in both the lattice spacing and the pion mass. The statistics needed to keep 
the errors in the flow scale fixed is found to increase approximately as the inverse square of the 
lattice spacing. Tree level improvement of the scales to and wo is found to be useful in most of the 
range of lattice spacings we explore. The scale uncertainty due to remaining lattice spacing effects 
is found to be about 3%. The ratio wo/\/to is Nf dependent and we find its continuum limit to be 
1.106 ± 0.007(stat) ± 0.005(syst) for mj^wo cs 0.3. 

PACS numbers: 


I. INTRODUCTION 

In any cutoff field theory it is easy to set the unit of mass in terms of the momentum cutoff. So, in lattice field 
theories the scale can be set by the inverse lattice spacing, 1/a. However, physically interesting questions require us 
to relate one measurable quantity to another, when both are computed to comparable precision in the theory. Using a 
physical scale to set the units of mass by eliminating the artificial choice of a is called setting the lattice scale. Doing 
this allows us to take the limit a —^ 0 in renormalizable theories without encountering artificial infinities. 

In principle, any mass scale can be chosen to define units, so the question of what to use for a mass scale is essentially 
one of convenience. An ideal scale should be easy to control numerically in the non-perturbative domain as well as 
be amenable to perturbative analysis. In recent years Wilson flow 0,i has emerged as a new and computationally 
cheap way of setting the lattice scale ii , since it seems to fulfill both criteria. 

However, Wilson flow scales, like Aj^, are theory scales. In order to determine them in “physical” (GeV) units, 
one needs two separate scale computations within the theory: one of the theory scale under question, the other of 
a measurable scale. Then by comparing the measurable scale to experiment, one can determine the theory scale in 
physical units. Glearly, in order to do this one needs to control two measurements. For Wilson flow this has been 
attempted in quenched QCD with 2 flavours of Wilson quarks 2+1 flavours of improved Wilson @,0 and 
improved staggered quarks Q and 2+1+1 flavours of improved staggered quarks Q. 

These computations have uncovered several systematics in the setting of the Wilson scale. In this paper we 
investigate in detail these systematics for two flavours of naive staggered quarks over a large range of lattice spacing 
and pion mass. We report on investigations of statistical uncertainties, as well as the dependence on all tunable 
parameters. We present an estimate of the Wilson flow scale in physical units. 

In the next section we outline the methods which we use. In Section nil we present a summary of the runs and 
statistics. A description of our results is given in Section HVl and a summary given in Section [Vj 
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II. METHODS AND DEFINITIONS 


Start with a gauge field configuration, i.e., the set of link matrices, {[/^(x)}, where x denotes a point in the 4-d 
Euclidean space-time lattice, and /i denotes one of the 4 directions. Wilson flow of this configuration is the evolution 
of these matrices in a fictitious “flow time” t, using the differential equation 


dUfj,{x,t) 

dt 


dS[U] 

dU^{x) 


Uf_,{x,t), where Uf_,{x,0) = Uf^{x), 


( 1 ) 


and the derivative is the usual Hermitean traceless matrix obtained by differentiating the scalar valued action func¬ 
tional 5'[f/] with respect to the link matrix Q. We use the convention that 


5[C/]=^Retr[l-C/(p)], 


( 2 ) 


where U{p) is the ordered product of link matrices around a plaquette, and the sum is over plaquettes. Clearly, the 
configuration with all [/ = 1 is a fixed point of the flow, and it can be shown that it is an attractive fixed point with 
a finite basin of attraction [I|. 

Following [l|, we define the scale by constructing the quantity 


£{t) = t^E{t), where 




( 3 ) 


where F^j^^ is a lattice approximation to the gluon field strength tensor and the bar denotes averaging over the lattice 
volume. The field strength tensor can be built either from the Wilson plaquette operator or through a 16-link clover 
operator. Some of our investigation of the systematics of Wilson flow involves comparing these two definitions. The 
scales which emerge from this are defined through the equations 


(^(^))lt=to(c) “ '"I ^ 


d{£{t)) 


dt 


= c. 


t=w^(c) 


( 4 ) 


The choice of c = 0.3 gives the quantities usually referred to as to and wq in the literature, a convention that we 
adopt. The modification, c=2l?> has also been suggested 0. The value c = 0.4 has been used in A weak coupling 
expansion [I| gives 


m) = 


(47r) 




If one uses to(c) to set the scale, then the expression above can be used to define a renormalized coupling 

IGtt^c 


9r = 


( 5 ) 


( 6 ) 


and it is clear that the choice of c is equivalent to a choice of the renormalization scheme. We report a study of this 
choice later in this paper. Note that the values of c used generally correspond to as = > 1. 

Tree-level improvement was performed by noting that the weak-coupling expansion in eq. m can be systematically 
corrected for lattice-spacing dependence through a computable piece 


m) = 


3 2c 
(47r)2^ V t 


where C 


t 


- 1 4- ^ C 2 


m—1 


t 


( 7 ) 


We use the coefficients presented in [1^. Later in this paper we show the effect of these corrections, and incorporate 
them in our measurements of the scale. 

We have also incorporated a finite volume correction due to the zero-mode of the gauge field m- Its effect is to 
scale 




1 - 


+ d(e-i/«') 




1 - 




( 8 ) 


where C = V^/L, L is the lattice extent, and d is a Jacobi Theta function. Except at our two smallest bare couplings 
the effect of the finite volume correction is comparable to, or smaller than, the statistical errors. 
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p 

ma 

Ns 

Lja 

Machine 

Traj 

(MD) 

Statistics 

To -1- T X A 

wola 

m-Ka 

5.2875 

0.1 

16 

V 

1 

400 -b 10 X 50 

0.6112 (4) 

0.790 (1) 


0.05 

16 

V 

1 

780 -b 10 X 50 

0.6354 (6) 

0.575 (1) 


0.025 

16 

V 

1 

200 -b 15 X 70 

0.6539 (1) 

0.415 (2) 


0.015 

16 

V 

1 

400 -b 10 X 50 

0.6608 (5) 

0.325 (2) 

5.4 

0.05 

16 

V 

2 

200 -b 20 X 75 

0.8418 (14) 

0.604 (2) 


0.025 

16 

V 

1 

400 -b 10 X 51 

0.9264 (21) 

0.443 (2) 


0.015 

24 

V 

2 

400 -b 10 X 50 

0.9600 (9) 

0.351 (1) 


0.01 

32 

G 

2 

200 -b 20 X 40 

0.9922 (7) 

0.292 (2) 

5.5 

0.05 

16 

V 

1 

200 -b 20 X 50 

1.1689 (40) 

0.613 (2) 


0.025 

24 

V 

1 

1680 -b 10 X 101 

1.2651 (18) 

0.446 (1) 


0.015 

28 

G 

2 

400 + 10 X 120 

1.3302 (13) 

0.353 (2) 


0.01 

32 

G 

2 

200 -b 20 X 40 

1.3771 (16) 

0.294 (2) 


0.005 

32 

BG 

1 

250 -b 10 X 50 

1.4254 (37) 

0.212 (1) 

5.6 

0.05 

24 

V 

1 

400 -b 10 X 55 

1.4850 (26) 

0.594 (2) 


0.025 

24 

V 

1 

1700 -b 10 X 48 

1.6007 (33) 

0.427 (2) 


0.015 

28 

G 

2 

400 + 10 X 120 

1.7087 (25) 

0.329 (2) 


0.01 

32 

G 

2 

200 -b 20 X 40 

1.7814 (36) 

0.272 (2) 


0.005 

32 

BG 

1 

300 -b 10 X 50 

1.8547 (71) 

0.198 (2) 


0.003 

32 

BG 

1 

600 -b 5 X 105 

1.8824 (32) 

0.151 (1) 

5.7 

0.025 

24 

V 

1 

530 -b 10 X 59 

1.9645 (48) 

0.395 (2) 


0.005 

32 

BG 

1 

370 -b 10 X 50 

2.1470 (73) 

0.177 (3) 


0.003 

32 

BG 

1 

300 -b 10 X 50 

2.2103 (162) 

0.139 (7) 


0.002 

32 

BG 

1 

480 -b 5 X 62 

2.3765 (67) 

— 


TABLE I: The data sets used in this paper. Runs were made at different bare couplings /3 and bare quark masses ma. 
Hypercubic Nf lattices were used, where Ns — L/a. N gauge configurations were collected for each run, after discarding 
an initial time Tq for thermalization, and collecting one configuration after every time T (both To and T are given in MD 
time units). The runs were performed on a vector machine (V), on a Blue Gene (BG) and GPUs (G). For completeness, our 
estimates of the tree-level improved flow scale wo/a and the pseudo-Goldstone pion mass m-na are also collected here, although 
they are discussed in detail in Section EYl 


III. RUNS 

We generated gauge field configurations with two flavours of naive staggered quarks over a wide range of bare 
couplings and bare quark masses. The bare parameters and statistics are given in Table ID Since the runs were 
performed on different machines we took the precaution of repeating several runs on multiple machines in order to 
cross check results. In these cases only the runs with the largest statistics are reported in the table above. 

A part of this range has been explored earlier, and pion masses have been reported [l^. We checked that at the 
common points our measurements of pion masses agree with those previously reported in the literature. All our 
analyses use the bootstrap technique to estimate expectation values and errors. Since we estimate pion correlation 
functions using a bootstrap, we reduce covariances between measurements at different distances by using independent 
bootstraps at each distance. Confidence intervals on the fitted parameters are estimated by a bootstrap over fits. 

The biggest challenge in estimating pion masses at small lattice spacings is in using lattice extents which are large 
enough to separate out the ground state from excitations. This is most acute for the Goldstone pseudoscalar mass at 
the smallest bare quark mass and lattice spacing, where our lattice size (mT^L < 4) was clearly inadequate. So we do 
not quote the pion mass from this lattice. Since we use naive staggered quarks, taste symmetry breaking remains a 
concern. We will report investigations of this elsewhere. 

One technical issue has to do with the integration of the flow equations. We tested both the Euler integrator and 
the fourth-order Runge Kutta (RK4) integrator. In Figure [1] we show the evolution of the plaquette under the flow 
when it is integrated using each of these methods for one fixed configuration. As shown, both integrators perform well 
even for dt = O.I. Note that with the Euler integrator and dt = 0.1 the first integration step has larger errors than 
the later steps. Such self-repair is seen also with a finer time step of dt = 0.01. With the Euler integrator one sees a 
failure of this self-repairing mechanism when dt = 0.2. The observation that the free configuration is an attractor of 



















4 



FIG. 1: The effect on the flow of changing integrator and dt. For dt < 0.2 the difference between integrators is negligible, 
especially with increasing t. For dt > 0.2 the integrator wanders off from the true solution. 


the map in eq. (HD serves to explain both self-repair and its failure. The fact that there is an attractor with a finite 
basin of attraction is the reason for self-repair, with the global errors being smaller than the 0{t^) predicted by a 
local analysis. Its failure occurs for sufficiently coarse dt, when the flow falls outside this basin of attraction. RK4 is 
generally more stable, but even so, its global error is smaller than local analysis would lead us to believe. We used 
RK4 with dt = 0.01, but checked the results statistically by changing dt by a factor of 4 either way. We found that 
the statistical uncertainty in the measurement of flow times is larger than any effect of the evolution. 

The statistical properties of the measurement of £(t) under evolution in flow time are also of interest. Since our 
measurements are separated by 10 or 20 MD trajectories, at t = 0 they are quite decorrelated. However, as the flow 
integrates information over successively larger volumes, one expects autocorrelations to grow with flow time. We 
quantify the autocorrelations in terms of the integrated autocorrelation time. Tint, which is defined in terms of an 
autocorrelation function of the measurements C{s) as 


T-int = 1 + 2 / dsC{s), (9) 

Jo 

where s is the separation between the measurements. In Figure [5] we show Tint for £{t) as a function of the flow 
time, t. In accordance with expectations, this shows an initial rapid increase. The observed plateau in Tint is due to 
insufficient statistics; clearly for the set with /3 = 5.6, the effective number of configurations decreases by a factor of 
around 20 when this plateau develops. Significantly more statistics would be needed to improve the measurement in 
this region, and decrease the estimate of the error in Tint elsewhere. 

If the scaling of autocorrelations is physical, i.e., has a sensible continuum limit, then a natural way to compare 
flow times for different simulations would be to scale them by iCq (or, equivalently, to)- HMC simulations with fixed 
trajectory lengths have Tint scaling as the square of correlation lengths [^. Since flow time also scales as the square 
of lengths, one should expect Tint — Figure [2] illustrates that this scaling is not present in the initial state, but 
develops fairly early during the flow and is a good first approximation to the observation. It would take significantly 
improved statistics to study the remaining deviations. The physics result is simple: as the lattice spacing decreases, 
the statistics required to keep a constant error on the flow scale increases (roughly) as the inverse square of the lattice 
spacing. 



5 



FIG. 2: The integrated autocorrelation time estimated for S(t) computed from the Wilson operator for several data sets. The 
autocorrelations increase very rapidly with flow time before running out of statistics. 



FIG. 3: The comparison of Wilson flow scales determined using the Wilson (plaquette) or the clover operator. The scale is 
varied by changing both gauge coupling and the bare quark mass. Up triangles denote v^q/o, down triangles are for wo/a. 
The unhlled symbols show the direct measurements, and filled symbols the tree-level improved values. The panel on the left 
collects all our measurements, that on the right shows only the four smallest lattice spacings, by zooming into the dotted box 
shown at the left. 


IV. RESULTS 

In Figure [3] we plot a flow scale obtained with the Wilson operator used for £{t) against the same scale obtained 
with the clover operator. If the two were equal, then the measurements would he on the diagonal line. It has been 
observed before that the clover improvement changes the flow scale i/to quite significantly, as we verify again. The 
data set for the scale wq is significantly closer to the diagonal. Both of these scales are improved significantly by a 
tree-level improvement, at least on coarser lattice spacings: both sets of measurements are moved significantly closer 
to the diagonal line. However, as shown in the zoom in Figure [3l the improvement is marginal for wq at the smallest 
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FIG. 4: The pion mass in physical units as a function of the bare quark mass. 


lattice spacings. This implies that any remaining finite lattice spacing corrections in wq are small. In view of this, we 
will use the tree-level improved value of wq to set the scale in the rest of the paper. We see that the range of lattice 
spacings we scan covers a factor of four from the coarsest to the finest. 

We apply this scale setting first to re-examine the pion mass measurement. Our measurements of in lattice 
units are given in Table ID We plot m,r in units of Wq in Figured] It is clear from the figure that the range of pion 
masses explored in this study covers a factor of four from the largest to the smallest. Given the rapid variation of 
a/wo and uIttWo with the bare coupling and the bare quark mass, it is useful to trade the bare parameters for these 
two. 

Since both the scales y/to and wg are physical, the ratio 7Z = wq/^/Iq is expected to tend to a good limit as 
the lattice spacing decreases. In Figure [S] we show the dependence of this ratio on the lattice spacing (given in 
units of the tree-level corrected value of wg)- At the smallest lattice spacing which we have examined {wg/a ~ 2.4), 
TZ. ~ I.100±0.003. For 2-I-1-I-I flavours of staggered quarks @ we deduce TZ ~ 1.21±0.01, where the error is estimated 
conservatively by neglecting covariance of the numerator and denominator. Since the statistical errors in TZ are small, 
the difference is significant. In a direct computation we checked that in the pure gauge theory, when wg/a ~ 2.4, the 
ratio TZ = 1.012 ± 0.005 (this is consistent with results presented in i). The ratio clearly depends on the number of 
flavours of quarks. 

TZ also depends on the lattice spacing and the quark mass, as shown in Figure O At fixed renormalized quark 
mass, rriT^wg ~ 0.3 we have tried a quadratic extrapolation to the continuum. Using the data points on the four finest 
lattices, the continuum extrapolated ratio is 1.101 ± 0.003. A fit using the quartic term gives the extrapolated value 
1.11 ± 0.01. If one uses only the three finest lattices, then the continuum extrapolation gives TZ = 1.108 ± 0.007. We 
put these observations together and quote a continuum extrapolated value 


= ^ = 1.106 ± 0.007(stat) ± 0.005(syst). 
v^o 

Following [^, we define a measure of the slope with respect to the lattice spacing as 

oa 7^(a = wo/1.75) 

= TZ{a = 0)-^ 


( 10 ) 


( 11 ) 


This is significantly larger than the results which can be reconstructed from values for other slopes quoted for Aly = 2 
clover improved Wilson fermions in [^ . At this time we are unable to comment on what combination of factors most 
influences this difference: the nature of the sea quarks, the value of m^, or technical issues in comparing slopes of 
slightly different quantities 0. 











FIG. 5: The ratio of the scales wo/y/to for all the bare couplings and quark masses studied. At the smallest lattice spacings 
the ratio from Wilson (up triangles) and clover (down triangles) operators, for both direct (unfilled symbols) and tree-level 
improved (filled symbols) scales tend to the same value. 




FIG. 6: (a) Variation of wo/y/to with the lattice spacing at fixed mnWo — 0.3. The systematics of the continuum extrapolation 
is discussed in the text, (b) Variation of wo/y/to with ttIttWo at the smallest bare coupling. Chiral logs M are not visible at 
this precision. 


It is known that wq is more strongly dependent on the quark mass than to 0 ] ■ A roughly linear dependence of both 
the scales with the renormalized quark mass has been observed before over a range of m,rWo similar to that explored 
here. Figure [6] shows this linear behaviour of the ratio wo/y/to. An extrapolation to the chiral limit as rri^WQ [1^ at 
our smallest bare coupling yields wo/y/to = 1.105±0.002 Using the Nf = 2 value for wq above. Defining an effective 
slope parameter 


QTTl _ 

— 


TZ{mTrWo 

TZ{mTrWo 


0.45) — TZ{mTrWo 
0.45) -I- TZ{mTrWo 


0.30) 
0.30) ’ 


( 12 ) 


our observations giv e ~ 2%. This is compatible with the change reported with two flavours of clover improved 
Wilson quarks in Q. 

Since the parameter c determines the value of the running coupling eq. one may use the RG-flow of the coupling 








FIG. 7: The ratio of the lattice spacings for two different sets of bare parameters depends on the renormalization scheme, i.e., 
the choice of c, through wo(c). 



P 



FIG. 8: The ratio extracted using the Lepage-Mackenzie scheme for computing . The first panel shows all 

our simulation results plotted as a function of the bare coupling j3. The next panel shows the results for fixed tUttWo — 0.3 and 
a quadratic continuum extrapolation. 


to examine the c-dependence of 'Wq{c). Define a measure of the change in Wo(c) through 


Wo(c = 0.4)/a — Wo(c = 0.2)/a 
2wo(c = 0.3)/a 


(a fixed). 


(13) 


On our finest lattice, we find Sc — O.I. The same measure with to gives about 0.2. The formal two-loop expression 
for the running of g/j in eq. (jH]) yields Sc — 0.3. Since the renormalized couplings obtained for these c are large, the 
two-loop beta function does not run the coupling reliably, so one should take the last number only as indicating that 
such large changes in scale are natural when changing c. 

It is more interesting to ask whether the ratio of lattice spacings at two different bare couplings and quark masses 
is independent of the choice of c, when each of these is given in units of Wq. Ideally, of course, such a ratio should not 
change with c. In Figure [7] we show that, in fact, there is some residual dependence on the parameter c. Take the 
case where this ratio is close to 1.175 for c = 0.3. The change in this ratio of lattice spacings for variation of Ac = 0.2 
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around c = 0.3 is 3% of the central value. While not ideal, this change is rather small. Presumably this uncertainty 
in the scale setting is due to remaining lattice spacing corrections. It would be interesting in the future to perform 
this comparison at smaller lattice spacings. 

Measurements of plaquettes can also be converted to a scale using the methods of [l^ . Since the scale setting by the 
flowed plaquette suffers from significant lattice spacing effects at flow times Wq, necessitating the various corrections 
which we have explored, it may be suspected that these effects could be larger at flow time t = 0. These are partly 
taken into account by corrections suggested in [T^. In Figure [5] we show the dimensionless ratio wq-^-ws obtained 
by a comparison of this scale with the flow scale wo/a. In the second panel of Figure [51 we show the ratio icoAms 
at fixed pion mass, uIt^wq ~ 0.3 as a function of the lattice spacing. One sees a strong, nearly quadratic, lattice 
spacing dependence, albeit with a slope smaller than A quadratic extrapolation to the continuum limit gives 
wqAj^ = 0.218 ± 0.001, where the error is statistical only. It is interesting to compare this indicative number to the 
value ioi Nf = 2 clover fermions. We take = 257 ± 26 MeV as quoted in [l^, and combine it with the value of 
wq reported with Nf = 2 clover fermions Q, to get wqAjjs = 0.23 ± 0.03. 


V. CONCLUSIONS 

We have reported on investigations of the Wilson flow scales ^/to and wq in QCD with two flavours of naive staggered 
quarks. Our investigations cover a wide range of lattice spacings (a factor of about 4) and pion masses (also a factor 
of about 4). We found that the scale wq has smaller lattice spacing artifacts than ^/io. One consequence of this is 
that tree-level improvement of the former has smaller effect than in the latter. In most of this paper we have used 
tree-level improved measurements of wg obtained from the clover operator as the object to set the scale by. 

We found an interesting approximate scaling of the autocorrelations of the basic measurement {S{t)). The integrated 
autocorrelation time increases with t before saturating. The scaling implies that keeping the error in measurements 
of wg fixed in the continuum limit may require the statistics to grow as the inverse square of the lattice spacing. 

We found that the ratio TZ = wg/^/tg = 1.100 ± 0.003 when wg/a ~ 2.4. A continuum extrapolation at fixed 
rriT^Wg ~ 0.3 gave TZ = 1.106± 0.007(stat) ± 0.005(syst). Comparison with results for the pure gauge theory, and with 
Nf = 2-I-1-I-1 reveals a dependence of TZ on Nf. The the compilation of Q also shows this trend for staggered 
quarks, but not for Wilson quarks. For Nf = 2 clover improved Wilson quarks, the value of TZ is different from our 
determination [^. 

The dependence of the scale wg(c) on c is large; this is natural since c enters linearly in the definition of gj^, which 
depends nearly logarithmically on the scale wg{c). In principle, this should not change the ratio of two lattice spacings. 
However, we found a mild (3%) dependence of the ratio of two lattice spacings on c. The effect is small enough that 
one suspects it is due to lattice spacing dependences which are not absorbed into the tree-level improvement of wg. 

By using our data sets to determine the scale via the Lepage-Mackenzie prescription [I^ we found that it has large 
lattice spacing corrections. However, with our data we tried a simple continuum extrapolation at fixed rriT^wg ~ 0.3, 
and found WgAjjg = 0.218 ± 0.001. If one then uses the ALPHA collaboration’s value Ajjg = 257 ± 26 MeV for 
Nf = 2, then one is led to the conclusion that for naive staggered quarks Wg = 0.17 ± 0.02 fm. This error is purely 
statistical, and dominated by the statistical error in the determination of Ajfs. 
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of the ILGTI in Mumbai, and with the Gray XK6 installation of the ILGTI in Kolkata. 
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